Theory of unconventional quantum Hall effect in strained graphene 
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We show through both theoretical arguments and numerical calculations that graphene discerns 
an unconventional sequence of quantized Hall conductivity, when subject to both magnetic fields 
(B) and strain. The latter produces time-reversal symmetric pseudo/axial magnetic fields (b). The 
single electron spectrum is composed of two inter-penetrating sets of Landau levels (LLs), located at 
±y2n|& ± B\, n = 0, 1, 2, • • •. For b > B, these two sets of LLs have opposite chiralities, resulting in 
oscillating Hall conductivity between and =F2e 2 //i in electron and hole doped system, respectively, 
as the chemical potential deviates from the neutrality point, but remains in its vicinity. The electron- 
electron interactions stabilizes various correlated ground states, e.g., spin-polarized, quantum spin- 
Hall insulators at and near the neutrality point, and possibly anomalous Hall insulating phase 
at incommensurate filling ~ B. Such broken symmetry ground states have similarities as well 
as significant differences from there counterparts in the absence of strain. For realistic strength of 
magnetic fields and interactions, we present scaling of interaction induced gap for various Hall states 
within the zeroth Landau level. 

PACS numbers: 71.10.Pm, 71. 10. Li, 05.30.Fk, 74.20.Rp 



Successful fabrication of two-dimensional electron gas 
(2DEG), e.g., Galium- Arsenide (GaAs) heterostructure, 
provided unique opportunity to observe a peculiar novel 
aspect of low-dimensional electronic systems, quantiza- 
tion of Hall conductivity {o~ xy ). At weaker magnetic 
fields (~ 1 T), even the samples with lower mobilities, 
discern quantized plateaus of a xy at various integer val- 
ues of e 2 /h. This phenomena is referred to as integer 
quantum Hall effect (IQHE) Rather more striking 
observation is the plateaus of that quantity at various, 
for example 1/3, fractions of e /h, in improved sam- 
ples, however at stronger fields (~ 10 T) [2|. Whereas 
the IQHE arises from free motion of fcrmions in mag- 
netic field Q, its fractional version necessarily requires 
strong electron-electron interactions to develop mobility 
gap within a partially filled Landau level [|[ . 

Integer quantization of a xy appears when the chem- 
ical potential (/z) lies within a mobility gap, filled by 
localized states, separated by two extended conducting 
edge modes, carrying the quantized Hall current As 
the magnetic field is reduced, more and more extended 
states get occupied and the total Hall current, which is 
the algebraic sum of the Hall current carried by each of 
them, encounters quantized increment, due to the iden- 
tical chirality of all the extended edge states, placed at 
well separated energies 0. 

Besides the GaAs heterostructure, the new generation 
two-dimensional electronic system, graphene, discerns a 
sequence of Hall plateaus at fillings v = ±4(n + \), sub- 
ject to relatively low fields Q, while additional plateaus, 
for example at v = 0,±1,±4, show up as the field is 
enhanced [8, 9]. Otherwise, all the LLs support the cur- 
rent carrying states with identical chirality, as in GaAs 
[lol ]. However, due to its mechanical flexibility under 
strain and the underlying hexagonal crystallographic ar- 
rangement, graphene may experience yet another mag- 



netic field, which may emerge from its bulging [11] . Such 
strain induced pseudo/axial magnetic field preserves the 
time reversal symmetry (TRS), and points in opposite di- 
rections at two inequivalcnt Dirac points, suitably chosen 
here at K = (1, l/\/3)(27r/a\/3) and —K [H[. Therefore, 
subject to strain as well as an external magnetic field, one 
can expose the gapless Dirac quasi-particles, near two 
Dirac points with different effective fields, \B ± b\, possi- 
bly pointing in opposite directions, respectively. Hence, 
an interplay of these two gauge fields, concomitantly an 
unconventional spectrum of the Hall conductivity can be 
realized in graphene. 

It is perhaps worth considering the Hall response of 
this system when B > b{^ 0) first, though such con- 
dition may not be achieved in laboratory at this mo- 
ment (l3l |. The spectrum of non- interacting Dirac quasi- 
particles is comprised of two inter-penetrating sets of LLs 
at well separated energies ±^j2n(B ± b) with degenera- 
cies (B ± &)/27r per unit area. A previous work predicts 
that quantized Hall conductivity then may appear at all 
the integer values of e 2 /h, however they occur at incom- 
mensurate fillings, due to distinct degeneracies of two sets 
of LLs [14j . The monotonic variation of the Hall conduc- 
tivity arises from the fact that all the LLs experience the 
orbital effect of the effective magnetic fields in the same 
direction, yielding the current carrying extended states 
with identical chirality. Hence upon increasing electron 
or hole doping, as the chemical potential sweeps through 
various LLs, the total Hall current adds up, leading to 
the announced spectrum of o~ xy . 

Rather more interesting situation arises when b > B. 
For B = 0, the pseudo Dirac LLs, spaced at ±y/2nb p~5l — 
[Lsl ]. and near two valleys have opposite chirality, hence- 
forth the TRS is preserved. As long as b > B > 0, tw o in- 
equivalent sets of LLs, now located at y2n(6 ± B), with 
respective degeneracies (6 ± B)/2n, per unit area, con- 
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tinue to enjoy opposite chirality (see left column of Fig. 
[I}. Consequently, as the chemical potential starts to de- 
viate from the charge neutrality point (CNP), the Hall 
conductivity is restricted within ±2e 2 //i (when more LLs 
localized near one valley is filled) and (equal number of 
LLs localized near two Dirac points are populated). The 
± sign corresponds to hole and electron doped systems, 
respectively, and note it is opposite to what one has in 
the absence of the strain-induced field, b. 

Even though the Hall conductivity stays bounded, as 
the chemical potential is enhanced, more and more cur- 
rent carrying edge states get filled. Therefore, in a 
two terminal conductance measurement, its longitudi- 
nal component G xx is expected to increase monotoni- 
cally. Every time the Hall conductivity suffers a quan- 
tized change, the two terminal longitudinal conductance 
G xx increases by 2e 2 jh (see right column of Fig. [TjlEoT]. 
However, the oscillatory sequence of o~ xy is strictly true, 
in the vicinity of the CNP. The spacing of the Dirac 
LL decreases with the LL index and the effective mag- 
netic field experienced by these sets of LLs arc different. 
Hence, far away from the CNP, LL crossing is unavoid- 
able, and the above mentioned spectrum of Hall conduc- 
tivity may cease. One can then in principle observe quan- 
tized plateaus of a xy at higher integer values > 3e 2 /h. 
When B <C b, the LL crossing occurs at large LL index 
(n). In that limit, assuming that the chemical potential is 
not too far from the Dirac points, one can safely neglect 
the LL crossing. Upon incorporating the Zecman split- 
ting, A z ~ B (T) K, additional Hall plateaus appears at 
a xy = ±e 2 /h, in between and ±2e 2 /h. 




FIG. 1: (Color online) Left: Spin degenerate interpenetrat- 
ing LLs of Hd[A, a]. Here we have shown the LLs for n — 
0, ±1, ±2 only. Two LLs have the degeneracies 2D± — (b±B) 
per unit area. Right: Schematic variation of Hall conduc- 
tances G X y (blue), and two terminal conductances G xx (red) 
as the chemical potential is tuned in a hole doped graphene. 
Here, E" — ^/2n(b + aB), with a — ±. In electron doped 
system, G„ remains the same, while G xy changes its sign. 

To compute the LL spectrum, we here construct an 
8-component Dirac spinor \P = ('J'-j-, v E , j,) T , where = 
[ut (K+q), »t (K+q), uU-K+q) , (-K+q)] , with a =t 
, i as electrons spin projection along the z— direction. 
The orbital effects of the real and pseudo magnetic fields 
can be captured by the Hamiltonian [T3 - [l6| 

H D [A, a}= I 2 ® «7o7i (ft ~ A ~ ^Ts^i) ■ (1) 
The magnetic fields are related to the vector poten- 



tials as B(b) — e^ijdiA(a)j. The gamma matrices are 
7o = h ® cr 3 , 7i = <r 3 <g) cr 2 , 12 = h® cri, 73 = u\ ® er 2 , 
7s = 02 ® c 2 [H} . The spectrum of Hp [A, a] is com- 
prised of two sets of interpenetrating LLs at energies 
±y / 2n(b + B) and ±y/2n(b — B), with respective de- 
generacies f2(6 ± B)/2n, shown in Fig. Q] [2l[. Here 
n = 0, 1, 2, ■ ■ ■ and f2 is the area of the strained graphene 
sample. With b > B states within the zeroth LL (ZLL) 
are localized on only one sub-lattice, say A for example, 
while they reside on complimentary sub-lattices near two 
Dirac points if B > b [Tj]- Otherwise, there are 2(6 ± B) 
states per unit area respectively near the Dirac points. 
The two-fold degeneracy of each LLs is due to fermion 
spin. The valley degeneracy for all the LLs at finite ener- 
gies are removed, as they are exposed to different effective 
magnetic fields. However, at precise zero energy LLs can 
be found near two Dirac points, guaranteed by an 'index 
theorem'. El Hi 




FIG. 2: (Color online) The energy spectrum (left) and wave- 
functions (WFs) (right) for a strained graphene in magnetic 
field when B < b. WFs localized on one edge, live on opposite 
side at two valleys, explicitly, A(B) and D (C) are localized 
on left (right) edge, therefore carrying opposite chirality. 

Let us firsts register the Hall response of the non- 
interacting system. At fi = 0, the ZLL is half-filled, 
containing total 4fi& states. A particle-hole symme- 
try, generated by I2 <8> 70 for example, anti-commutc 
with Hr>[A,d\\ of the spectrum then guarantees exact 



zero Hall conductivity, o x 



0. Upon doping, when 



the chemical potential is placed between the ZLL and 
yb — B, o~ xy remains at zero. It can be seen easily after 
subscribing the Stf eda formula [25| for the Hall conduc- 
tivity <j xy = (§7j)^, in the natural units e = c = 1. A is 
the electronic density in the bulk below the chemical po- 
tential. Derivative with respect to B is taken at fixed /x, 
measured from the half-filled band. In order to place the 
chemical potential such that < /x < \fb — B, one needs 
to fill N = 2flb (independent of B) states from the CNP, 
yielding a zero Hall conductivity. On the other hand, if 
\/b— B < n < y/b + B, 5N = -25 B and hence a xy = 
—2. The factor 2 counts the spin degeneracy of the LL . 
Upon further doping when y/b + B < fj, < \/2(b — B), 
the Hall conductivity comes back to zero. Hence, with 
odd (even) number (modulo 2 due to spin) of LLs below 
the chemical potential, one gets a xy = —2(0), as long as 
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there is no LL crossing. The underlying reason for such 
oscillatory nature of the Hall conductivity is the follow- 
ing: two sets of LLs near ±K experience total effective 
magnetic fields (b±B) > 0, which point in opposite direc- 
tions and hence the current carrying extended states of 
these two LLs have opposite chirality. When odd number 
(modulo 2) of LLs above the CNP are filled, a xy = —2, 
since there is imbalance in the occupation of the LLs 
near two Dirac points. With an even number (modulo 2) 
of filled LL a xy = 0, as the Hall currents from the LLs 
near +K and —K exactly cancels each other. With hole 
dopping, the Hall conductivity oscillates between and 
+2 e 2 /h. 



from all the LLs, including the ZLL, \Jn(b ± B) — > 
\J n(b ± B) ± A z . As the Zeeman gap scales as A z ~ B 
K (<C J2a=± a \/ n (b + °~B)), if B is measured in Tesla, 
it does not cause LL crossing near the CNP. Otherwise, 
the Hall conductivity remains pinned to zero, when fi = 
due to the particle- hole symmetry, generated by o\ <g> 70, 
for example. It remains so even when the chemical poten- 
tial lies in between the Zeeman shifted ZLL at A z and 
\/b — B — A 2 , because 5N = 2bil (independent of B). 
Otherwise for v& + &B + oA z < fi < \/b + oB — oA z , 
a xy = - (^ s -) e 2 /h, where a = ±1. Therefore, the 
Zeeman splitting introduces additional Hall plateau at 
=Fe 2 //i, in electron and hole doped systems, respectively. 
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FIG. 3: (Color online) Two possible interaction driven split- 
ting of the zeroth Landau level and spectrum of the Hall con- 
ductivity, a xy — fe 2 /h, at various chemical doping. ± corre- 
sponds to the ZLL states localized near the valley at ±K. 

The chiral nature of the edge modes in the presence 
of strain and magnetic field can also be seen in a finite 
honeycomb lattice, with zigzag edge. The orbital effect 
of the real magnetic field can be captured by attach- 
ing a Pcicrls phase, £y — > tij exp (27ri-^ f? A ■ all) to the 
nearest-neighbor (NN) hopping amplitudes. The strain 
induced axial gauge fields can be incorporated in the sys- 
tem by local modulation of the NN hopping amplitudes. 
We here modify the hopping only along one of the three 
bonds, oriented orthogonal to the zigzag edge [U [22|. 
With such simple deformation, even though one ends up 
with slightly inhomogeneous Fermi velocities and thus 
LL engergies, it is still sufficient to capture the pecu- 
liarities of the edge modes arising from the preserved 
TRS. One can see from Fig. [5] that the chiralities of 
two states localized near one zigzag edge of the system 
arc opposite, if b > B, same when B = 0, but b =^ 0.[2lj] 
However, the LLs near two Dirac points appears at dif- 
ferent energies. Therefore, as one changes the chemical 
potential the Hall conductivity keeps oscillating between 
±e 2 /h and (since we consider spinless fermion here). If 
on the other hand, the real magnetic gets stronger, i.e., 
B > b, the edge modes near two Dirac points share iden- 
tical chiralities. [2LJ Therefore, when the chemical poten- 
tial varies, o~ xy changes monotonically. [l4j An interesting 
possibility is B — fe, keeping only one of the Dirac points 
exposed to finite magnetic field, yielding plateaus of o~ xy 
at v = 2(n + i). The other one remains semi-metallic, 
giving finite o~ xy and o~ xx simultaneously. [2li [26j 

The Zeeman splitting (A z ) lifts the spin degeneracy 



Due to the peculiar entanglement of the ZLL with the 
sub-lattice degrees of freedom (all the states are local- 
ized on one sub-lattice), a strain induced charge den- 
sity wave order (not spontaneously generated) always 
persists within the ZLL. Nevertheless, this configuration 
is a natural ground state for the residual NN Coulomb 
repulsion. Since two valleys are populated by distinct 
number of states, a 'valley polarized' quantum anoma- 
lous Hall state cannot develop at CNP. It may however 
be realized at incommensurate filling, v oc B about the 
neutrality point. Upon occupying two valleys with op- 
posite spin projections, a spin Hall state, corresponds 
to circulating currents among the sites of same sub- 
lattices of honeycomb lattice, can be realized. It also 
carries a finite ferromagnetic moment, oc B (difference 
of LL degeneracies). Since the ZLL is localized on one 
sub-lattice, the ferromagnct (FM) order is tied with an 
anti-ferromagnet (AF) order. The spin Hall order pa- 
rameter reads as A QS h = (^jcr (8 ijijo^fn-). ^. l28j In 
the absence of the gauge fields, though it requires finite 
strength of the second neighbor interaction (V2) to trig- 
ger such ordering [2^], even an infinitesimal V2 can de- 
velop a finite spin Hall order, when the kinetic energy of 
the Dirac quasi-particlcs is quenched by the axial fields 
[H-[l8j and to the leading order A QSH oc V 2 & The 
spin polarized (SP) state carries finite FM (Afm) and 
AF {A AF = (*Jct <g) 7o*<t)) orders simultaneously [l]|- 
Though the Zeeman effect yields Afm 7^ 0, the onsitc- 
Hubbard interaction (U) .which is possibly the strongest 
interaction in graphene [3l|, prefers an AF ordering, in 
the vicinity of the CNP. A similar spin polarized state 
can also be realized even when B = [171 ]. However, in 
Ref. [TtJ SP state was identified as pure ferromagnet, 
which appears to be limited in light of above discussion. 
The AF order parameter anti-commutes with Hd[A,o\. 
Hence, apart from splitting the ZLL, it optimally low- 
ers the energy of the filled Dirac sea, by p ushing all the 
LLs at finite energies further down, ±y/2n(b ± B) — ¥ 
±^2n{b±B)+A 2 AF . The spin-polarized gap within the 
ZLL is A S p = A AF + A FM - U. 

For small Zeeman coupling one can write the ground 
state energy per unit area at half filling in the presence 
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of AF (N) and spin Hall (C) orders as 

-. N 2 C 2 
E[N,C] = - + - + E o[ N,C], 



(2) 



where g a { c ) ~ U (V2). E [N, C] is the ground state energy 
per unit area of the effective single particle Hamiltonian 

Hhf = H D [A, a] - (N ■ a) <gi 70 - (c ■ &\ ® 17172- (3) 

With small Zceman coupling spin-anisotropy can be ne- 
glected and we set N(C) = (N(C), 0,0), for simplicity. 
Minimization of E[N, C], with respect to N and C, yields 
two coupled gap equations, which for N > C reads as 
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for i,j = 1,2, where g = (g a ,g c ),X = (N,C),£ = (b,B) 
and e„, a = [2n(6 + aB) + (N + aC) 2 ] 1 ' 2 . The first term 
in the gap equations are ultra-violet divergent and inde- 
pendent of Xi. However, cut-off independence of the gap, 
demands g a = g c for both AF and spin Hall order to be si- 
multaneously finite. Since, in graphene, U > V2, possibly 
a spin polarized state (TV = Aaf 7^ 0,C = Aqsh = 0) 
is formed at the CNP. Even though, with b > B, there 
exists a series of a xy = plateau, only the one near \i = 
bears an AF order, while the rest arises from lack of 'val- 
ley reflection symmetry'. Placing the chemical potential 
close to the first excited state at ±A SP , a spin Hall order 
develops additional incompressible Hall state, leading to 
a xy — ±e 2 /h ( see left column of Fig. 3). If on the other 
hand, V2 > U, yielding C > N, the splitting of the ZLL 
gets reversed ( Fig. 3, right column). 




FIG. 4: (Color online) Dimensionless activation gap (A/A) for 
Cxy = 0, ±e 2 //i (black, red) Hall states when b = 300 T and 
T < B < 50 T, for 5 = 0, 0.03, 0.07, 0.14, 0.3, 0.7 (from top 
to bottom), assuming A ~ 1/2. 5A, the ultra-violet cutoff over 
which the dispersion is approximately linear. g c is the zero 
field criticality. Lower x-axis denotes B/A 2 (dimensionless). 

Minimizing the ground state energy, one can find the 
gap equation for a xy = plateau near the CNP. For fixed 
axial magnetic field (b = 300T), the interaction induced 
gap at the neutrality point increases linearly with the 
real magnetic field (0 T < B < 50 T) (See Fig. 0}, since 
with increasing B LLs at lower energies (— e n! +) is pushed 



further down and also their degeneracies (~ b + B) in- 
creases. Otherwise the gap is insensitive to exact nature 
of the order parameter. Scaling of the activation gap for 
a xy = ±e 2 /h state within the ZLL is qualitatively same 



as (however smaller than) that for o~ xy = 0. Such hierar- 
chy comes solely from the ZLL, as exactly half of the ZLL 
contributes to the gap for a x 



state, while a fewer 



states within the ZLL contribute to the <r xy = ±1 state. 
f2ll ]. Otherwise, activation gap for both the Hall state, 
scales sub-linearly with the strength of the interaction, 
S = (gA)^ 1 — (g*A) -1 , where g* is zero field criticality 
for insulation. If either of the magnetic fields becomes 
inhomogeneous, the LLs at finite energies get distorted, 
however the ZLL stays unaffected. Therefore, interaction 
induced gap generation take place even when the field is 
non-uniform, whereas the gaps then closely follow pro- 
file of the magnetic fields. [24[ With weak inhomogeneous 
fields, the quantization of o~ xy is expected to survive. 

In the absence of the axial field or even when B > b, 
the states within the ZLL, localized near two valleys live 
on complimentary sub- lattices [3, [32j . Therefore, a con- 
ventional AF order develops by filling up states on two 
sublattices with opposite spin projections, if the onsite 
repulsion is strong enough. However, the staggered spin 
moments on two sub-lattices are of different magnitudes. 
Therefore, one may argue such a correlated phase as fer- 
rimagnetic as well [l4j . On the other hand, when b > B, 
the AF state is formed by filling states in ZLL, with a 
particular spin projection, living on one sub-lattice. 

In experiment 11|, the axial field is localized only in 
certain region of the sample, however it is uniform. Par- 
ticles circling the axial field inducing bumps can pick up 
an axial Aharonov-Bohm phase (ABP), only if it travel 
through the strained region. Otherwise, it is identical for 
the trajectories in opposite directions, whereas the ABP 
due to the real magnetic field are of opposite sign for the 
trajectories in opposite directions [333 • Consequently, the 
trajectories with opposite circulation, acquires different 
effective ABP, namely, sum and difference of it due to 
two fields. Therefore, in Hall conductivity measurement, 
the terminals need to be attached to the regions with at 
least finite strain, since the axial gauge potential is pro- 
portional to strain, though b can be zero. In molecular 
graphene^^, and recently fabricated strained graphene 
on Ru substrate [35( , the axial field can possibly be real- 
ized in the entire sample, thereby may easily discern the 
peculiar Hall conductivity, we propose here. 
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Supplementary material for "Theory of unconventional quantum Hall effect in 

strained graphene" 



I. DIAGONALIZATION OF DIRAC 
HAMILTONIAN IN REAL AND AXIAL 
MAGNETIC FIELDS 

We here present one possible prescription to diagonal- 
ize the Dirac Hamiltonian in the presence of real and 
strain induced pseudo/axial magnetic fields. Recall, the 
low energy Hamiltonian in graphene, subject to strain 
and magnetic field reads as [l[ 

H D [A, a]= h<8> «7o7i (& —Ai — «73750i) . (5) 

The axial vector potential (a*) is a member of general 
SU(2) time reversal symmetric gauge potential 

a-i = 73a? + 75a- + ij3l5ai 5 - (6) 



However, the first two terms break the translational sym- 
metry, generated by I tr = I 2 <£> 27375 • If the deformation 
of the graphene flake is smooth enough one can safely 
neglect their contributions and only keep the af 5 compo- 
nent. Upon setting a| = = 0, we set af 5 = a i; for 
notational simplicity. Thereafter, the Dirac Hamiltonian 
is devoid of any valley mixing. One can then exchange 
the second and the third 2x2 blocks of Hd[A, a] and 
cast it in a block diagonal form, say H+ © H- , where 

H± = ±I 2 ®cT 1 (qiTai-A 1 ) + l2®cT2(q2Ta2-A 2 ). (7) 

H± capture the effect of both the gauge potentials on the 
Dirac fermions in the vicinity of ±K points, respectively. 
However, both H± are unitarily equivalent to the cele- 
brated Dirac Hamiltonian, Hp = «7o7i(<7i — A e /), subject 
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FIG. 5: (Color online) Graphene with zigzag edges. The edge 
(along x direction) is effectively infinite and the width in y 
direction is finite. In the strained case, we suppose the hop- 
ping terms along y direction has a gradient which ranges from 
t — 1 — W/2 to t = 1 + W/2 from one side to the other. 




FIG. 6: (Color online) The low energy spectrum (left) and the 
local density of states (LDOS) (right) for isotropic graphene 
sheet which has infinite length in x direction and finite width 
of TV = 300 unit cells in y direction. The zero energy edge 
states exist between two Dirac points and contribute peaks in 
the LDOS. 



ia2 and A e J = a, + A, , 
i<J2®l2 and AS- = ai—Ai 



to effective magnetic fields, b±B respectively. Explicitly, 
H+ = U\H D U U with Ui=I 2 < 
while H- = U ] 2 HU 2 -, with U 2 

0,0. After expressing the original Hamiltonian Ho [A, a] 
as direct sum of two inequivalent copies of the standard 
Dirac Hamiltonian subject to different effective magnetic 
fields, one can immediately diagonalize it, yielding the 
announced spectrum of inter-penetrating sets of Landau 
levels (LLs) 



II. NUMERICAL DIAGONALIZATION OF 
NON-INTERACTING HAMILTONIAN IN 
STRAINED FINITE SIZE GRAPHENE, SUBJECT 
TO MAGNETIC FIELDS 

We consider a single layer graphene with infinite length 
in x direction and finite width of TV unit cells in y direc- 
tion as shown in Fig. [5] The tight-binding Hamiltonian 
for spinless fcrmions (for simplicity) with only nearest- 
neighbor (NN) hopping reads as 



ilbj + H.c), 



(8) 




FIG. 7: (Color online) The energy spectrum, LDOS and wave- 
function in a real magnetic field. In this case, the Landau 
levels are clearly developed. The Fermi energy crosses the 
first Landau level at four points A, B, C, and D, in which A, 
B and C, D belong to two different valleys. 



where <Zj and bj are the fermionic annihilation operators 
on sites of two inter-penetrating triangular sub-lattices. 

stands for the summation over three NN sites. In 
pristine graphene, the hopping amplitudes along three 
directions are equal. The above Hamiltonian can be 
represented by a set of Harpers equation, due to the 
translational symmetry along the x-direction. One can 
then easily obtain the energy momentum relation numeri- 
cally. The energy spectrum and the local density of states 
(LDOS) for isotropic zigzag edge is shown in Fig. [BJ The 
Dirac points are located at k x = 2tt/3 and k x = 47t/3. In 
the vicinity of the Dirac points the dispersion is linear. 
The edge states at E = arise due to charge accumula- 
tion on two different sub-lattices, causing peaks in LDOS, 
which possibly can be detected in STM measurement. 

The orbital effect of the real magnetic field can be cap- 
tured by introducing Peierls phase in the NN hopping 
amplitudes, Uj — > tije 2m ^ ij . 4>ij is related to the mag- 
netic vector potential as 



— / A ■ dl. 

he 



(9) 



In the presence of magnetic field, the linear dispersions 
around the Dirac points quenche into sets of Landau lev- 
els at well separated energies, as shown in Fig. [7] It is 
evident from this figure that the wave functions living 
on the same edge share identical chirality. For example, 
wave-functions on left side of each valleys, A and C, are 
located on left zigzag edge, whereas those on the right 
side of each valleys, B and D, are located on right zigzag 
edge. 

One can introduce the axial magnetic field (b) by lo- 
cally modifying the hopping amplitudes. One particular 
way of such deformation is the following. We here modify 
the hopping only along one of the three bonds, oriented 
orthogonal to the zigzag edge. It smoothly varies from 
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FIG. 8: (Color online) The energy spectrum (left) and the 
wave functions (right) in a strain field. The hopping, perpen- 
dicular to zigzag edge, ranges from t - W/2 to t + W/2 with 
W = 1.0 from one edge to the other. Here we set t = 1. Wave- 
functions localized on one edge, live on opposite side at two 
valley, explicitly, A(B) and D (C) are localized on left (right) 
edge. Henceforth, the edge states near two valleys have differ- 
ent chirality. In the presence of real magnetic fields, situation 
is reversed, i.e., wave-functions on left (A,C) (right (B,D)) 
side of each valleys are located on left (right) edge. 




FIG. 9: (Color online) The energy spectrum when B ~ b. 
The Landau levels are only developed in one valley and the 
Dirac cone is recovered at the other one. 



t - -f to t 



where W < t, from one edge of the sys- 



tem to the other, separated by 300 unit cells. With such 
simple deformation, one ends up with slightly inhomogc- 
neous Fermi velocities and thus LL energies. Neverthe- 
less, we can immediately see from Fig. [5]that chirality of 
the edge modes are opposite at each edge in the presence 
of axial/pseudo magnetic field as we discussed in the 
manuscript. 

Upon tuning the relative strength of the fields, one 
can achieve an interesting limit, b ~ B. Only one of the 
Dirac points is then exposed to finite magnetic field, 2B 
or 2b. Therefore, the energy spectrum near that Dirac 
point is composed of Landau levels. On the other hand, 
the system remains semi-metallic in the vicinity of the 
other Dirac point [f|. This scenario is depicted in Fig. |H1 



III. DIRAC LANDAU LEVELS WITH 
SYMMETRY BREAKING TERMS 

Let us now compute the LL spectrum of the emergent 
single particle Hamiltonian, containing various symmetry 
breaking terms. For our purpose, we keep both the anti- 
ferromagnet (./V) and the quantum spin Hall (C) orders 
besides the free Dirac Hamiltonian with real and axial 
vector potentials. The resulting Hamiltonian reads as 



Hhf = h <E> «7o7» (?« - M - hzlsai) - [N ■ a) <g> 70 
- (p ■ a) <g> hi 72 + X z (<r 3 <8 h) , (10) 

if one wishes to keep the Zeeman coupling (A z ) as well. 
However, in comparison to the broken symmetry order 
parameters, the Zeeman term is much smaller and one 
can neglect it. Then all spin projections are equivalent 
and we choose N = (N, 0, 0) and C = (C, 0, 0) for the 
ease of calculation. Even in the presence of symmetry 
breaking terms two inequivalent valleys are decoupled 
and one can cast Hhf is block diagonal form H + © H_ , 
where 



H± = ±I 2 ® cri (<?i T ai - M ) + I 2 



oiiki T« 2 - A 2 ) 
(11) 



after exchanging the second and third 2x2 blocks. Using 
the unitary rotations U± and U2, defined in sec. I, one 
can cast H± as 

H ± ^i l0ll (q l Ta l -A l )- l0 (N±C). (12) 

Therefore the spectrum of Hhf is composed of two sets of 



inequivalent LLs, at energies ±y 2n(b ± B) + (N ± C) , 
for n > 1 and ±|iV± C\ when n = 0. 

IV. GAP EQUATIONS 

The ground state energy per unit area at half-filling in 
the presence of two symmetry breaking order parameters 
reads as 



E\N,C\ = 



AT2 
45a 



^9c 



E [N,C], 



(13) 



as shown in the main text. Eq[N, C] is the ground state 
energy of the effective single particle Hamiltonian Hhf- 
At Half-filling, the ground state of Hhf has all the states 
at negative energies filled, while those at positive energies 
at completely empty, yielding 

E [N,C] = ((b + B)e n ,+ + (b-B)e n ^ 

77 n>l ^ ' 



— Ub + B){N + C) + (b-B)\N-C\ , 

(14) 
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where e„ i(T = \J 2n(b + aB) + (N + aC) , for a = ±1. In 
the main text of the paper, we argued that unless g a = 
g c , the anti-ferromagnct and the spin Hall order cannot 
coexist. For g a > g c , i.e when the on-site repulsion (U) is 
stronger than the next neighbor one (V2), at half- filling 
the ground state configuration is with while C = 

0. Minimizing the ground state energy, E[N, 0], with 
respect to N, then yields the gap equation at half filling 



1 



1 



47r 3/2 
b 

2-KTTl ' 



s 3/2 



2s 



er=±l 



e 2s(b+aB) _ l 



(15) 



where g = g a . The first term on the R. H. S. counts 
the contribution of all the filled LLs with n > 1, and it is 
ultra-violet divergent, while the second term on the R. H. 
S. arises only from the ZLL. However, the divergence can 
be regulated by defining the zero field critical interaction 
for ordering as 



1 



2A 



ds 

772 : 



(16) 



after ir/g — > 1/g, so that the physical gap remains cut- 
off (A) independent. The cut-off (A) defines the range of 
energy, over which the dispersion is approximately linear. 
Defining the strength of the interaction as 



<y = _i L 

Ag Ag* ' 



(17) 



where (5 measures the deviation of the interaction from 
the zero field critical interaction for insulation, one can 
cast the gap equation as 



x 1 



S_ 

m 



1 



1 



1 



f{x). (18) 



Here we introduced new variables, x = ^tp-, and q = 
and the function f(x) reads as 

m = 1 



1 



ds 
^371 



2se~ s l x 
e 2s - 1 



- 1 



(19) 

Similarly, one can arrive at the gap equation for the 
a xy = ±e 2 /h Hall state, which on the other hand reads 



as 



2.S 



4tt 3 / 

b + B 
4to7t 



e 2s(b+aB) _ I 

(20) 



where g = g c , assuming N > C ov g a > g c . Here the 
contribution from the filled LLs (first term on R. H.S.) is 
identical to that in Eq. (Tf"5")) . but the ZLL contribution 
is different from that in Eq. (|15p , since only fewer states 
contributes to the ground state energy when the chemical 
potential is tuned close to the spin polarized gap (N). 
After performing similar exercise, one can finds the gap 
equation to be 



m 0r V 



(21) 



For 8 = (g = </*) and b — 300 T, solutions of the 
gap equations are at x = xq = 5.0537 and 17.9226 for 
B = T, respectively. For B = 50 T, the solutions are 
at x = 5.741 and 15.378.0 For a finite 8 > (g < g»), 
the solutions can only be found at x > xq. Hence, for 
realistic strength of two fields, one can safely expand f(x) 
for large argument of x, leaving us with simple algebraic 
equations for the interaction induced activation gap for 
&xy = 0,±e 2 //i Hall state, read as 



x 1 



1 

71 



s_ 

m 



-^ + 0(x 
\ x 



-3/2n = 



) =0,(22) 



8_ 

m 



( 1 + ^f) (^«+^ + ^- 3/2 ))-». 



(23) 



respectively. The coefficients are u = 2.0652 and v = 
0.924. The numerical solutions of the gaps are shown in 
the main part of the paper, for a set of realistic strengths 
of interactions (S > 0) in graphene. 
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